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ABSTRACT 

Shallow-water numerical simulations show that the atmospheric circulation of the close-in extrasolar 
giant planet (EGP) HD 209458 b is characterized by moving circumpolar vortices and few bands /jets (in 
contrast with ~ 10 bands/jets and absence of polar vortices on cloud-top Jupiter and Saturn). The large 
spatial scales of moving circulation structures on HD 209458 b may generate detectable variability of the 
planet's atmospheric signatures. In this Letter, we generalize these results to other close-in EGPs, by 
noting that shallow- water dynamics is essentially specified by the values of the Rossby {Ro) and Burger 
{Bu) dimensionless numbers. The range of likely values of Ro (~ 10~^-10) and (~ 1-200) for the 
atmospheric flow of known close-in EGPs indicates that their circulation should be qualitatively similar 
to that of HD 209458 b. This results mostly from the slow rotation of these tidally-synchronized planets. 

Subject headings: planetary systems - planets and satellites: general - stars: atmospheres - turbulence 



L INTRODUCTION 

The focus of extrasolar planet research has broadened 
to now include the characterization of their physical prop- 
erties, as shown by the recent sodium detection in the 
atmosphere of HD 209458 b (Charbonneau et al. 2002). 
Atmospheric circulation is expected to play a key role in 
determining a number of observational characteristics of 
EGPs, including their albedo and transmission spectrum 
(see, e.g., Seager & Sasselov 1998, 2000; Sudarsky et al. 
2000; Brown 2001). This is especially true for close-in 
EGPs, which arc thought to be tidally-locked to their par- 
ent star and irradiated on one side only: circulation will be 
essential in redistributing heat from the day to the night 
side on these planets, thus determining to a large extent 
how they will appear to the distant observer (Cho et al. 
2002a,b; Showman & GuiUot 2002). 

Recently, we have presented a set of detailed shallow- 
water numerical simulations of the atmospheric flow on 
HD 209458 b (Cho et al. 2002a,b), currently the only 
EGP with known mass and radius from the transit light 
curves and radial velocity measurements (Charbonneau et 
al. 2000; Henry et al. 2000; Mazeh et al. 2000; Jha et 
al. 2000; Brown et al. 2001). These simulations suggest 
that, contrary to the simple day/night (hot/cold) picture, 
the circulation on this planet is characterized by two mov- 
ing circumpolar vortices and a small number of latitudi- 
nal bands /jets. The vortices act as dynamically distinct 
thermal spots whose motion around the poles generates 
variability as seen by an observer interested in quantities 
integrated over the planetary disk (or circumference). 

It is possible to determine the general features of the 
circulation pattern expected within the framework of 
shallow-water dynamics by specifying the two dimension- 



less numbers - Rossby [Ro) and Burger {Bu) - for the 
atmospheric flow. In this Letter, we estimate a range of 
likely Ro and Bu values for known close-in EGPs and con- 
clude that their atmospheric circulation pattern should be 
qualitatively similar to that of HD 209458 b. In §2, we re- 
call how the atmospheric flow pattern can be characterized 
by the knowledge of Ro and B^,- In §3, we describe the 
sample of close-in EGPs selected for our study and how we 
estimate likely values for various global planetary param- 
eters entering into the deflnition of Ro and Finally, 
our results and conclusions are presented in §4. 

2. TURBULENT SHALLOW- WATER DYNAMICS 

Shallow- Water equations describe the motion of a thin, 
homogeneous layer of hydrostatically-balanced, inviscid 
fluid with a free surface, in motion around a rotating 
planet (Pedlosky 1987, Holton 1992). The fluid is subject 
to gravitational and Coriolis forces and obeys the following 
equations 

V • Vv ^ -gVh - /k X V, (1) 



dt 
dh 

at 



+ V • V /i 



-W- V, 



(2) 



where v is the horizontal velocity, h is the thickness of 
the modeled layer, / — 2f2sin(/3 is the Coriolis parame- 
ter, n is the rotation rate of the planet, tp is the latitude, 
g is the gravitational acceleration and k is the unit vec- 
tor normal to the surface of the planet. In dimensionless 
form, shallow-water equations become functions only of 
the Rossby [Ro) and Burger (i?„) numbers: 



Ro — 



U 



(3) 
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where U, L and H are characteristic velocity, length and 
layer thickness scales, respectively; L^i is the Rossby de- 
formation radius. Note that |/| ~ f2 at mid-latitudes and 
that the planetary radius, Rp, is the relevant length scale 
when discussing the large-scale atmospheric circulation. 
The Rossby number measures the importance of rotation 
on the flow, while the Burger number measures the strati- 
fication of the atmosphere via the Brunt- Vaisala frequency 
(Holton 1992). 

The atmospheric structure in bands of gaseous giant 
planets in our Solar System is well described as emerg- 
ing from freely-evolving shallow-water turbulence on the 
sphere (Cho & Polvani 1996b). Turbulence in a thin at- 
mospheric layer is quasi-2D in nature. Contrary to the 
forward turbulent energy cascade observed in 3D geome- 
try, 2D turbulence is characterized by an inverse energy 
cascade (transfer from small to large scales) and a forward 
cascade of enstrophy"^ down to small scales, where it is 
dissipated by viscous processes. Qualitatively, the inverse 
cascade is associated with the growth of vortices through 
continuous mergers. 

Turbulence in a thin atmospheric layer is also strongly 
constrained by the combined eS'ect of spherical geometry 
and rotation (the effect"). While the force balance 
is everywhere the same along a latitude circle, it changes 
with latitude because of the dependence of the Coriolis 
term with (p. This anisotropy, as measured by the param- 
eter /3 = 2Q cos (fi/Rp (the latitudinal gradient of /), is 
strongest at the equator, where /3 is maximum.. While 
fluid motions arc free to grow to the largest available scale 
in the longitudinal direction, their growth is limited in 
the latitudinal direction by the characteristic Rhines scale, 
Lfj = 7r-\/2[///3 (Rhines 1975). This anisotropic growth 
is a likely origin of the banded structure on gaseous giant 
planets in our Solar System; the number of bands expected 
for a given planet is roughly A^band ~ i^Rp/Lfj. 

Cho & Polvani (1996a) presented an extensive numerical 
study of freely-evolving shallow-water turbulence. They 
explored the entire parameter space of the equations, as 
determined by the two dimensionless numbers Ro and Bu- 
They showed that the anisotropy due to the /3— effect on a 
rotating sphere is necessary but not sufficient to produce 
a long-lasting banded structure. The Rossby deformation 
radius, Ld, must also be ^ i?p/3 for the banded structure 
to be stable. This small value of Lu, which acts as a lim- 
iting scale for vortex interactions, prevents the formation 
(via successive mergers) of large-scale structures such as 
circumpolar vortices. 

Cho et al. (2002b) presented a generalization of the re- 
sults of Cho & Polvani in the case when the planet is sub- 
ject to day-side hemispheric heating, as expected for close- 
in EGPs. Day-side heating was prescribed in the adiabatic 
limit by forcing the fluid to be permanently thicker on that 

^The flow cnstrophy is defined as , where ^ = k ■ V X v is the flow vorticity. 

*We hinit the eccentricity of planets in our sample to e < 0.1, by analogy with the small eccentricities of Solar System giants to which 
"zero-eccentricity" shallow- water models have been applied with success (Cho & Polvani 1996b; Cho et al. 2002b). 
^http : //exoplanets . org/almanacf rame . html 
®http : //www . obspm . f r/encycl/ encycl . html 

'^The periastron distance, where tidal forces are the strongest, is also roughly the circular radius expected for the orbit after complete 
circularization. 



side, while keeping the average thickness constant. Exten- 
sive exploration of the parameter space of the forced model 
showed that previous shallow- water dynamics results were 
recovered even in the presence of this extra forcing (i.e. L/j 
and L n remain the relevant scales determining the number 
of bands and the formation of polar vortices). 

These results can be recast in terms of the values of 
Ro and i?„ for the atmospheric flow. By setting L ~ Rp, 
I/I ~ n and /3 ~ ^/Rp (mid-latitudes), we see that the 
number of bands/jets expected is A^band "^1/ \/{2Ro) and 
that the presence of circumpolar vortices is expected for 
Bu ^ 1/9. Thus, if the values of Ro and B^ for other 
close-in EGPs can be estimated, one can get an idea of 
the type of large-scale atmospheric circulation expected 
on these planets in the stable, radiative region. In our 
estimates of Ro and B^, we will set J7 = C/, which is the 
global velocity scale of the atmospheric flow and is only 
known for the Solar System giants (Table 1), 

An important parameter entering the definition of both 
Ro and B^ is the planetary rotation rate, fi. While the 
value of Q. is generally unknown for EGPs, a number 
of close-in EGPs have the advantage of being probably 
tidally-synchronized to their parent star, so that their ro- 
tation rate has effectively been measured via the orbital 
period (O = Oorb for circular orbits). As we show be- 
low, the knowledge of for this sample of close-in EGPs 
restricts the range of possible values for Ro and B^ to 
a small enotigh region of the parameter space that their 
atmospheric circulation pattern can be inferred. 

3. SAMPLE OF CLOSE-IN EXTRASOLAR GIANT PLANETS 

Since; tidal synchronization occurs faster than orbital 
circularization, it is possible that some close-in EGPs 
with substantial eccentricities are nonetheless (pseudo- 
)synchronized (i.e. synchronized at the periastron orbital 
frequency). The shallow- water results described in §2 were 
established only in the limit of negligible eccentricity, how- 
ever. Hence, we must restrict our sample to planets with 
small eccentricities.^ Table 1 lists all the EGPs selected 
for our study, plus the four Solar System giants. Parame- 
ters for the EGPs were collected from the extrasolar planet 
almanac^ and encyclopedia.^ 

Low-eccentricity EGPs were divided into two groups, 
based on their orbital distance to the parent star. In the 
first group, EGPs with semi-major axes a < 0.066 AU 
are most likely tidally-synchronized since all known EGPs 
with such small values of a are also circularized. The tidal 
synchronization status of the more distant planets (in the 
second group) is less clear because several eccentric EGPs 
with distances of closest approach^ as small as 0.05 AU are 
also known. It is thus not clear why some EGPs with pe- 
riastron distances larger than this value would be tidally- 
circularized while others would not. We note, however, 
that for values of the tidal parameter Q not too different 
from that of Jupiter (~ 10^), EGPs in this second group 
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are also expected to be synchronized. We will assume it is 
indeed the case in our calculations. 

Radial velocity surveys only measure MpSini, which is 
a lower limit to the planet's mass, Mp, given the unknown 
orbital inclination, i. For randomly oriented systems, the 
distribution of cos i is uniform. We adopt the value of Mp 
corresponding to sini = 0.5 for our fiducial estimate of 
Ro and and we allow sini to vary from 0.1 to 1 when 
estimating the range of likely values for Ro and 

For a given mass, Mp, the radius of an isolated planet 
is estimated from the mass- radius relation for sub-stellar 
objects of Chabrier & Barafi^e (2000), supplemented at the 
low mass end by a constant density law that empirically 
fits values for the Solar System giants. To account for the 
slower cooling under strong stellar irradiation, we also al- 
low the radius to be up to 50% larger than the value for 
an isolated planet, in agreement with published cooling 
EGP models (Burrows et al. 2000; Guillot & Showman 
2002). Our results depend only weakly on the planetary 
radius, as long as Rp ^ Rjup (as expected for all masses 
of interest). The gravitational acceleration is derived as 
g = GMp/Rp, where G is the gravitational constant. 

For the mean layer thickness, H, we adopt the atmo- 
spheric pressure scale- height, -ffatm = "T^-Tktm/.?, where TZ 
is the perfect gas constant. The global radiative equilib- 
rium temperature of the planet is Tatm = T*(-R*/2a)^/^(l — 
ylfc)^/^, which is a function of the parent star luminosity 
(Li, oc T^R^), the planet's semi-major axis a, and Bond 
albedo Af,. We adopt Ai, = 0.5 for all our numerical esti- 
mates; our results only weakly depend on the value of Ab 
unless it approaches unity. The stellar luminosity is de- 
rived from the mass through the simple mass-luminosity 
relation = {M^/Mq^-^Lq. 

The last two parameters needed to determine Rg and i?„ 
are the planetary rotation rate O and the global kinetic en- 
ergy scale tj . We assume that $7 = Oorb (as determined 
by radial velocity surveys) in all cases. We allow U to 
vary from 50 m s""'^, the smallest observed value for giant 
planets in the Solar System (Jupiter), to 1000 m s~^, a 
rather large value for which the typical wind speeds in the 
atmosphere of hot, close-in EGPs approaches the sound 
speed. A value U = 400 m s^^ is adopted for our fiducial 
estimate of Ro and B^. 

4. RESULTS 

Estimated values for Rg and _B„ arc given in Table 1 for 
Solar System giants and close-in EGPs. The values listed 



for close-in EGPs correspond to the range of min./max. 
values found given the various assumptions detailed in §3. 
Fiducial estimates are also reported in figure 1 , where solid 
dots correspond to group 1 EGPs (safe tidal synchroniza- 
tion assumption) and open circles to group 2 EGPs (tidal 
synchronization assumption less safe). HD 209458 b is 
indicated as a star. 

It is clear from figure 1 that close-in EGPs, as a group, 
occupy a different region of the Ro-Bu parameter space 
than Solar System giants. In particular, they systemat- 
ically have a Burger number B^ > 1/9 (even when ac- 
counting for the large range of allowed values; Table 1), 
which indicates that the presence of circumpolar vortices 
is expected in the radiative region of close-in EGPs within 
the framework of shallow-water dynamics. The larger val- 
ues of Ro also indicate that generally few bands/jets are 
expected on these planets (the uncertainty on U strongly 
affects this number; see Table 1), thus allowing the forma- 
tion of larger "great spots" (which could also contribute to 
the variability: Clio et al. 2002b). The near alignment of 
all the points representing close-in EGPs in figure 1 shows 
that the dominant parameter determining their position in 
this diagram is their rotation rate {Ro oc Bu oc f7~^). 
The small values of Ro and Bu for Solar System giants 
reflect their relatively fast rotation rates. 

Although we argued in favor of variable atmospheric sig- 
natures for close-in EGPs, it is important to note that 
models do not yet quantitatively predict how much vari- 
ability is expected. InChqetal. (2002a, b), we emphasized 
that the combination of U (unknown) and the amplitude 
of day-night heating (parametrized in adiabatic simula- 
tions) determines the contrast of the thermal spots asso- 
ciated with circumpolar vortices. In the future, diabatic 
shallow-water models will allow a self-consistent determi- 
nation of the day-night forcing. More sophisticated mod- 
els, combined with detailed radiative transfer and chem- 
istry descriptions (Seager & Sasselov 1998; 2000; Seagcr 
et al. 2000), will allow us to make r[uantitative predic- 
tions regarding the level of variability expected for various 
atmospheric signatures. 
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Table 1 

GLOBAL PLANETARY PARAMETERS 



Planet 


Mi, 


^orb 


a 


e 


Mp 






9 


n 




H 


u 


Ro 






(Mq) 


(days) 


(AU) 




("Jup) 




(nT) 


(m s- = ) 






(m) 


(m s-1) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


(8) 


(9) 




(10) 


(11) 


(12) 


(13) 


Neptune 


1.0 


60,189.0 


30.05 


0.0113 


0.054 


2 


5 X 107 


11 


9.75 X 10~^ 


3 


X lO* 


300 


1.2 X 10^1 


5-6 X 10~^ 


Uranus 


1.0 


30,685.4 


19.20 


0.0457 


0.046 


2 


6 X lo"^ 


9 


(-)l X 10-* 


3.r 


i X 10* 


300 


1.2 X 10^1 


4-8 X 10^2 


Saturn 


1.0 


10,759.2 


9.58 


0.0565 


0.30 


6 


X lo'' 


9 


1.6 X 10~* 


4 


X 10* 


300 


3.2 X 10-^ 


3.9 X 10~3 


Jupiter 


1.0 


4,332.6 


5.20 


0.0489 


1.00 


7. 


.1 X lo'' 


23 


1.8 X 10"'' 


2 


X 10* 


50 


4.0 X 10~3 


2.8 X 10~3 


Group 2 


rho CrB b 


0.95 


39.6450 


0.2300 


0.0280 


>1.100 








1.8 X 10~8 








0.24-7.4 


68-180 


HD 195019 b 


1.02 


18.3000 


0.1400 


0.0500 


>3.430 








4.0 X 10~^ 








0.10-4 


30.8-69.6 


Gl 86 b 


0.79 


15.7800 


0.1100 


0.0460 


>4.000 








4.6 X 10~8 








0.09-3.4 


13.2-46.4 


55 Cnc b 


0.95 


14.6530 


0.1150 


0.0200 


>0.840 








5.0 X 10~^ 








0.084-2.8 


13.2-39.2 


HD 130322 b 


0.79 


10.7240 


0.0880 


0.0480 


>1.080 








6.8 X 10~S 








0.062-2 


6.8-18 



Group 1 



HD 168746 b 





.92 


6.4090 





.0660 





.0000 


>0.240 




1 


.1 


X 


10-^ 




0.036-1 


.92 


3.2-22.4 


HD 49674 b 


1 


.00 


4.9400 





.0600 





-0000 


>0.120 




1 


.5 


X 


10~^ 




0.030-1 


88 


2-4-23-6 


Ups And b 


1 


.30 


4.6170 





.0590 





-0340 


>0.710 




1 


.6 


X 


10~^ 




0.026-0 


96 


2-44-8-0 


51 Peg b 


0. 


.95 


4.2300 





.0512 





-0130 


>0.440 




1 


.6 


X 


10~^ 




0.024-1 


.04 


1.64-7.6 


HD 209458 b 


1 


.05 


3.5247 





.0450 





.0000 


0.690 


10* 


8 2 


.1 


X 


10~^ 


2 X 10^ 


0.024-0 


.48 


3.6 


HD 75289 b 


1 


.05 


3.5100 





.0460 





.0540 


>0.420 




2 


,1 


X 


10-^ 




0.024-0 


.88 


1.32-6.0 


BD -10316 b 


1 


10 


3.4870 





-0460 





.0000 


>0.480 




2 


,1 


X 


10-^ 




0.020-0 


84 


1.36-6.0 


Tau Boo b 


1 


30 


3.3120 





-0500 





-0000 


>4.090 




2 


2 




10~^ 




0-019-0 


72 


1.36-4.96 


HD 179949 b 


1 


,24 


3.0930 





.0450 





-0500 


>0.840 




2 


,4 


X 


10-5 




0.018-0 


.62 


1.20-3.56 


HD 187123 b 


1 


.06 


3.0900 





.0420 





.0300 


>0.520 




2 


.4 


X 


10~^ 




0.019-0 


.72 


1.08-4.4 


HD 46375 b 


1, 


.00 


3.0240 





.0410 





.0000 


>0.249 




2 


.4 


X 


10-5 




0.017-0. 


.90 


1.0-6.4 


HD 83443 b 





.79 


2.9853 





.0380 





.0800 


>0.350 




2 


.4 


X 


10-5 




0.017-0 


.80 


0.80-6.4 



NOTES: (1) In order of decreasing orbital period (2) Stellar mass (3) Orbital period (4) Semi-major axis (5) Eccentricity (6) 
Planet mass (7) Planet radius (8) Surface gravity (9) Rotation rate (10) Atmospheric scale-height (11) Global atmospheric 
velocity scale (12) Rossby number (13) Burger number 
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Fig. 1. — Location of Solar System and close-in extrasolar giant planets in the Rossby-Burger space. The assumption of tidal synchronization 
for extrasolar planets represented by solid circles is the safest (group 1; see Table 1). HD 209458 b is indicated by a star. A representative 
range of possible values around the fiducial estimates for extrasolar giant planets (each individual solid or open circle) is shown as an errorbar 
(see Table 1 for details). Formation of circumpolar vortices is expected in the region to the right of the vertical dotted line (Bu > 1/9). A 
larger number of bands is expected for smaller values of Ro (see text for details) . 



